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ABSTRACT 


This report presents the results of comparing four different frequency estimation schemes in the 
presence of high dynamics at low carrier-to-noise ratios. 

The comparison is based on measured data from a hardware demonstration. The tested algorithms 
include a digital phase-locked loop, a cross-product automatic frequency tracking loop, an extended 
Kalman filter, and finally, an FFT-aided cross-product frequency tracking loop. The tracking 
algorithms are compared on their frequency error performance and their ability to maintain lock 
during severe maneuvers at various carrier-to-noise ratios. 

The measured results are shown to agree with simulation results carried out and reported previously. 
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Chapter 1 


Summary 


1.1 Executive Summary 

Results from this task have established that several algorithms for tracking high-dynamic Global 
Positioning System (GPS) signals meet the tracking requirements of the GPS Range Application 
Joint Program Office (RAJPO) [1]. These algorithms include several frequency trackers developed 
at the Jet Propulsion Laboratory (JPL) as well as the frequency tracker used by Interstate Elec- 
tronics Co. (IEC) in RAJPO’s Translator Processing System (TPS). All the algorithms presented 
here were evaluated via software simulations, and key algorithms were evaluated with hardware in 
real-time using a breadboard version of NASA’s Deep Space Network (DSN) Advanced Telemetry 
Receiver (ARX). 

The algorithms were tested in an environment which simulates the high dynamics and low carrier 
signal-to-noise ratios (CNR) that are typical of critical mission stages encountered in RAJPO 
applications. The highest dynamic stress is an acceleration ramp of 50 g in 0.5 s, or jerk (derivative 
of acceleration) of 100 g/s. This is an upper bound on the performance of some modern agile 
missiles. The CNR scenario requires tracking well below 30 dB-Hz, representing reception of GPS 
signals with antenna gain of -10 dB and a receiver front-end noise figure of 3.5 dB [2]. This low 
CNR is often forced by the small physical size of the host vehicle and by possible nulls, or low-gain 
zones, in the receiving antenna. Under such high dynamics and low CNR conditions, typical GPS 
receivers are incapable of tracking the Doppler frequency of the carrier. The algorithms presented 
here, as well as the algorithm developed and demonstrated by JPL in a preceding task [2], are 
capable of maintaining tracking under these conditions with a frequency error of a few Hertz. Once 
the high-dynamic maneuver is completed, carrier tracking by the receiver’s phase-locked loop (PLL) 
can be resumed. 

A related issue is the initial acquisition of GPS signals under high dynamics. Acquisition is usually 
measured in terms of Time To First Fix (TTFF), which is the time from receiver warm-up and 
signal presence until a reliable position solution is available. In some missions a TTFF of less 
than 5 seconds is required, either to satisfy range safety needs or to accommodate vehicles with 
very short mission durations. The TPS accomplishes a low TTFF by utilizing a parallel search 
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in the code-lag domain. We present results here which show that even lower TTFFs are feasible 
via parallel searches in both the code-lag and frequency domains. The latter utilizes Fast Fourier 
Transform (FFT) processing. 


1.2 Objectives and Approach 


The first objective of this task is to expand the high-dynamics tracking techniques developed under a 
preceding task [2] and to investigate and validate other related tracking and acquisition techniques. 
A second and equally important objective is to transfer the resulting technology to application in 
the industrial sector through RAJPO’s contractor, IEC. 

To accomplish these objectives, several algorithms were developed and documented. Key algorithms 
were: 


1. Maximum Likelihood Estimator (MLE) — an FFT-based method that is an extension of the 
method demonstrated in the preceding task; 

2. Extended Kalman Filter (EKF) — a phase and frequency tracker that utilizes near-optimal 
recursive estimation techniques; 

3. Automatic Frequency Control (AFC) Loops — in particular we consider three types of loops: 
a Cross-Product AFC (CPAFC), an Overlapping Discrete Fourier Transform AFC (ODAFC), 
and a Frequency Extended Kalman Filter (FEKF); 

4. FFT-CPAFC — an algorithm developed by IEC and used in the TPS; 

5. Adaptive Least Squares (ALS) — algorithms that apply optimum estimation techniques to 
the tracking of bi-phase shift keying (BPSK) modulation signals, e.g., when GPS data-aiding 
is not available; and 

6. Fast acquisition algorithms using FFTs. 


These algorithms were evaluated via analysis and software simulation, and performance comparisons 
were derived. Some of these algorithms were also evaluated on a breadboard version of NASA’s 
ARX, hence demonstrating their real-time practicality. 

Finally, progress in this work was reported to RAJPO, and through RAJPO to IEC, so as to ac- 
complish technology transfer. The TPS, under development for RAJPO by IEC, incorporates a 
FFT-CPAFC loop and meets RAJPO’s present acquisition and tracking requirements. Other tech- 
niques presented here are applicable to future equipment that may require even higher dynamics, 
lower CNRs, or shorter TTFF. 
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1.3 Results 


Algorithm performance was verified via analysis, simulations, and real-time hardware experiment 
where possible. Results of the three evaluation approaches are in close agreement, and exceed 
RAJPO’s requirements for the TPS. 

The tracking algorithms were extensively evaluated using two trajectories. The first trajectory has 
maximum dynamics of 100-g/s jerk, implemented as a 50-g change in 0.5 s, and repeated twice, 
while the second trajectory, used in testing the FFT-CPAFC (IEC’s algorithm), has maximum 
dynamics of 70-g/s jerk, implemented as 70-g change in 1 s. Both trajectories represent more 
severe dynamics than RAJPO’s specification of 50 g/s for 1 s. Results were compared in terms of 
root-mean-square (rms) frequency error and probability of loss-of-lock. 

Results are presented here in terms of CNR at the output of the code correlators which, when 
code lock exists, are modeled as 2-ms integrators. Since RAJPO’s documents specify the CNR at 
the output of the GPS receiving antenna we assume in this report that a 3-dB loss occurs in the 
system ahead of the code correlator output, i.e.,in the quantizers, bandpass filters, code correlator, 
etc. Hence, RAJPO’s typical CNR specifications of 31 dB-Hz and 38 dB-Hz appear here as 28 dB- 
Hz and 35 dB-Hz, respectively. This assumption must be included in all comparisons to external 
specifications for consistency. 

Figure 1.1 presents simulated rms frequency error as a function of CNR for the DPLL, MLE, FFT- 
CPAFC, EKF, ODAFC, FEKF, and CPAFC algorithms with 100-g/s jerk. RAJPO’s specification 
for maximum allowed rms error, 51.2 Hz (32 ft/s) at 28 and 35 dB-Hz and 50-g acceleration, is 
clearly met by all the algorithms evaluated. The conversion factor from ft/s to Hz is 1.6 Hz/(ft/s), 
based on the Doppler effect at GPS LI frequency. All the algorithms simulated present virtually no 
loss-of-lock at 28 dB-Hz under the specified dynamics. One must note that RAJPO’s specification 
is for 50-g dynamics with unspecified dynamics before and after the maneuver, while our results 
are for the particular trajectory defined later. 

Acquisition algorithms were evaluated for various dynamics. The critical dynamics parameter here 
is the maximum (unknown) velocity, which corresponds to the size of the frequency domain area 
that must be searched. Results show that RAJPO’s specifications, TTFF of 5.3 s at CNR of 28 
dB-Hz and 3.9 s at 35 dB-Hz for dynamics of 6150 ft/s, are met by the two-dimensional search 
algorithm, and that the portion of TTFF allocated to channel acquisition can be reduced to less 
than 0.5 s. The part of TTFF that is allocated to settling of the navigation Kalman filter is not 
affected. 


1.4 Conclusions and Recommendations 


Results presented here show that the high-dynamics tracking algorithm implemented in the TPS 
meets all of RAJPO’s requirements. Other algorithms presented here provide improved perfor- 
mance in terms of both tracking and acquisition capabilities. These should be considered for 
implementation for future suitable applications. 
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Figure 1.1: RMS frequency errors for a 100-g/s trajectory 





Chapter 2 


Purpose and Scope 


2.1 Purpose 


The purpose of this final report is to document the GPS high-dynamic tracking concepts that have 
been developed under the task plan, to evaluate the performance of these methods, and to present 
performance comparisons between the different algorithms, including the FFT-CPAFC loop used in 
the TPS procurred by RAJPO. This report fulfills deliverable item (e) of JPL’s Task Plan RE-182, 
Amendment 452A, “High Dynamic GPS Range Instrumentation Receiver” [3], sponsored by the 
U.S. Air Force System Command, Armament Division. 


2.2 Scope 


This report documents the concepts, testing, results, and conclusions per the task plan. Chapter 
3 presents a historical perspective of this task and a definition of the high-dynamic environment 
model used in the remainder of this report for both simulations and measurement. Chapter 4 
briefly describes the various high-dynamic tracking algorithms developed at JPL and IEC. (Detailed 
algorithm descriptions are included in the references.) Chapter 5 describes the experimental setup 
used at JPL to evaluate real-time performance. Chapter 6 presents results of algorithm evaluation 
both in software simulations and using the experimental setup. Chapter 7 concludes this report by 
highlighting some of the significant results of this task. 
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Chapter 3 


Introduction 


3.1 History 

JPL has pioneered the development of high-dynamic parameter estimation techniques for GPS 
receivers and related applications. This effort began in 1983 when JPL received a contract from 
RAJPO to validate a proposed concept for tracking high-dynamic vehicles without the use of 
inertial aiding. The problem arises because high dynamics introduce correspondingly high Doppler 
frequency shifts on the GPS RF carrier signals. When these Doppler shifts are large enough, the 
receiver’s PLL (actually a Costas loop) cannot maintain lock, and carrier tracking is lost. Since 
carrier tracking provides the navigation solution with accurate measurements of range rate, the 
quality of the range and velocity estimate is degraded. In addition, the loss of carrier tracking 
prevents recovery of the 50-Hz data modulation. The JPL concept was to apply an MLE approach 
to carrier frequency tracking, using an FFT processor. This approach allows for recovery of range 
and range rate estimates even in the presence of very high dynamics. Range can also be recovered 
from the code delay locked loop. 

In 1984-85, JPL developed a breadboard single-channel simulation system and verified performance 
of this concept with actual hardware. The design was mostly digital so that a proposed receiver 
architecture was suitable to miniaturization using VLSI technology. Breadboard tests showed that 
both pseudorange and range rate could be tracked simultaneously in the presence of 50-g to 100-g 
accelerations at CNRs as low as 28 dB-Hz, without inertial aiding. These results were also confirmed 
via analyses and extensive simulations. Results of that work are documented in [2] and [4]. 


3.2 Outline of Completed Research 

Subsequent JPL activity carried out under this task plan focused in two areas: tracking of high- 
dynamic GPS signals in the presence of data-wipe and tracking of BPSK-modulated signals. Data- 
wipe consists of the use of a secondary channel to remove the 50-Hz data modulation from the 
GPS signal. This can be performed in RAJPO’s environment when a translator is used (Figure 
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3.2.1). In this case, the high-dynamic vehicle carries a frequency translator that receives the GPS 
signals in L-band, translates them to a different frequency band, and then transmits them to a 
ground station. The GPS signal is also received from the GPS satellites directly at the ground 
system and the 50-Hz data is recovered using a separate GPS receiver. The recovered 50-Hz data 
is then applied to the re-transmitted signal (from the translator) with appropriate delay to “wipe” 
the data modulation. This approach, if properly implemented, results in the reception of a Doppler 
unmodulated tone, rather than a harder-to-track BPSK-modulated sine wave. The second area of 
JPL’s work was improved tracking of BPSK-modulated signals using differential techniques. 

Several frequency trackers were evaluated for the data- wipe case. These included extensions of the 
MLE concepts proven in the preceding task, a novel application of the EKF, ramifications of the 
CPAFC loop, and a DPLL. The fine tuning and performance of these frequency estimators were 
shown to depend on the signal dynamics, and were related to the maximum allowable observation 
time and highest-order derivative of the frequency process. Software simulations showed that both 
the MLE and the EKF can track trajectories that include 100-g/s jerk pulses with estimation errors 
of only a few Hertz at CNRs as low as 24 dB-Hz. In addition, the MLE maintained lock with 90 
percent probability over this severe trajectory even at 23 dB-Hz. 

After simulations of these algorithms were completed [5], several of these algorithms were demon- 
strated using a modified version of NASA’s Advanced Receiver breadboard as a testbed. As the 
Advanced Receiver breadboard was not originally designed for this purpose, several hardware and 
software modifications were required. These modifications are described in Chapter 5. The PLL, 
CPAFC, EKF and a hybrid algorithm developed by IEC, denoted an FFT-CPAFC algorithm, were 
implemented and tested on the breadboard, confirming the analyses and simulations. 

In the second area, tracking of BPSK-modulated carriers, work was based on various ALS al- 
gorithms. These algorithms use optimal estimation techniques and are rather computationally 
intensive, but provide superior performance. 


3.3 Definition Of High-Dynamic Trajectory 


A fundamental question in developing techniques for the tracking of high-dynamic vehicles is “What 
are high dynamics?” Usually, one interprets this as high values for velocity and its derivatives. 
In the preceding task [2], two high-dynamic trajectories were used: one that simulated circular 
motion, or turns, the other simulating linear acceleration. In the first trajectory the velocity and 
its derivatives are sine waves, corresponding to motion in a circle with a period of 6 to 8 s, and 
with radial acceleration of 50 g. This trajectory exibits an infinite number of derivatives and is 
hard to track with an algorithm that assumes a polynomial model for the velocity. The second 
trajectory assumed that the acceleration is constant throughout, except for a step change of 50 g. 
Both trajectories represented hard-to-track environments. 

In this task, the performance of the various estimators was evaluated on the basis of their ability 
to track a common trajectory. The trajectory chosen for both the simulation and the hardware 
demonstration experiments is derived from RAJPO’s specification [1] and uses step jerk (derivate 
of acceleration) as the most stressing dynamic. The FFT-CPAFC, used in the TPS, was evaluated 
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GPS SATELLITES 



Figure 3.2.1: Tracking scenario using GPS frequency translators (figure courtesy IEC) 


with a trajectory that included a single 1-s, 70-g/s jerk step, for a total acceleration change of 70 
g. This trajectory is shown in Figure 3.3.1. The other algorithms were evaluated with a higher- 
dynamic trajectory that included positive and negative jerk pulses of 0.5-s duration and 100-g/s 
magnitude, separated by 2 s of constant acceleration, as shown in Figure 3.3.2. The velocity 
trajectory can be scaled to an equivalent Doppler-frequency trajectory using the conversion factor 
of 5.25 Hz/(m/s), for LI signals. 

The dynamics affecting an individual tracking channel depend on the relative geometry of the GPS 
system. Each channel “sees” the component of dynamics that is along the line of sight between 
the receiver and the satellite assigned to this channel. Since this component is always less than or 
equal to the overall dynamics, the stress on each channel is not as severe as indicated by the overall 
dynamics, i.e., in a 50-g turn some channels may “see” no acceleration while others respond to an 
acceleration of some magnitude less than or equal to 50 g. Nevertheless, in this report we assume 
the worst case: the dynamics are along the line-of-sight from the receiver to the GPS satellite. 

The common input to the algorithms simulates the output of the code correlators sampled every 
2 ms. These samples represent in-phase and quadrature sinusoids derived from the frequency 
trajectory and observed in the presence of noise. All algorithms were evaluated on the basis of 
their response to these trajectories. Each 2-ms sample consisted of a complex number defined with 
real and imaginary components, and as follows: 


Xi = Acos(0 t ) + n Xt (3.1) 


Vi = Asin(0,) + n yt (3.2) 

Therefore, (a:**, 2/* ) is the complex sample (i.e., the in-phase and quadrature samples), 9{ is the ith 
phase sample, T s is the sampling time (2 ms), and n Xi and n y% are additive white Gaussian noise 
samples each with variance a 2 . The CNR is defined by: 

cnr = 2 Sr 3 dB ~ Hz (3 - 3) 

As mentioned before, this is the CNR at the output of the code correlator, or integrator. When 
comparing results in this report to other GPS data, one must translate the CNR at the correlator 
output to CNR at the RF input with the appropriate system loss. 
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Figure 3.3.1: Simulated 70-g/s trajectory 
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3.3.2: Simulated 100-g/s trajectoi 
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Chapter 4 


Algorithm Descriptions 


In this chapter, various tracking and acquisition algorithms are presented. Descriptions of most of 
the algorithms have already been published in the literature in one form or another, so no attempt is 
made here to reproduce the analysis. Instead, a brief description is given and the reader is directed 
to the appropriate publications for more detailed information. Algorithms are presented in three 
different categories. First, we present frequency tracking algorithms for data-wipe environments, 

i.e, MLE, EKF, DPLL, CPAFC, ODAFC, FEKF, and FFT-CPAFC. Then, algorithms for frequency 
tracking that apply ALS algorithms are presented. Finally, algorithms for fast initial acquisition of 
GPS signals are discussed. 


4.1 Maximum Likelihood Estimator (MLE) 

The structure and performance of the MLE of frequency in the presence of dynamics and additive 
noise is detailed in [5]. Here we present a brief description of the basic operations carried out by 
this estimator and summarize its performance by means of estimated rms estimation errors and 
estimated loss-of-lock probabilities. 

The MLE bases its estimate of frequency (and its time derivatives) on a vector of N consecutive 
in-phase and quadrature samples. The estimator structure is based on the maximization of the 
likelihood function, derived from the conditional joint density of the observed samples, conditioned 
on the relevant signal parameters. It is shown in [5] that even for severe dynamics characterized by 
100-g/s jerk pulses, it is sufficient to carry out the maximization of the likelihood function over the 
two-dimensional frequency-frequency rate plane, without incurring large estimation errors near the 
operating threshold. Thus, to first order, the effects of jerk (the rate of frequency rate) can often 
be ignored. However, the simultaneous estimation of both frequency and frequency rate is crucial, 
even though only the frequency estimates are of interest. 

The maximum likelihood estimator can be implemented by evaluating the likelihood function over 
a grid of points in the frequency-frequency rate plane, with spacing fine enough to resolve the likeli- 
hood function near its peak. This implementation can be reduced to a sequence of FFTs performed 
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on a suitably weighted input sequence, where the complex weights depend on the frequency- rate 
coordinate. Rough estimates of frequency and frequency rate are simply the coordinates of the 
peak, which can be refined by means of interpolation algorithms as needed. 

The simulation results are summarized in [5]. The rms estimation errors exhibit a well-defined 
threshold, below which a rapid increase in the estimation error occurs. This is also the region 
where loss-of-lock begins to be a problem (loss-of-lock is a condition in which the frequency estimates 
become independent of the true frequency). Here we define the loss-of-lock threshold as that CNR 
at which the estimator loses lock ten percent of the time. For N = 80 (N is the number of 
samples used in the FFT before zero-padding), simulation results were also obtained using the true 
value of jerk, in order to bound the possible improvement in the frequency estimates and loss-of- 
lock performance if the jerk was estimated as well. Apparently, no improvement is possible near 
threshold. Above threshold, however, rms frequency estimation errors could be reduced by the 
simultaneous estimation of jerk, acceleration, and velocity components. With the given trajectory, 
a CNR threshold of 23 dB-Hz was achieved with N = 50, resulting in a 7-Hz rms error. 


4.2 Extended Kalman Filter (EKF) 


Kalman filters have been shown to be the optimum estimators when the received signals are lin- 
ear functions of the unknown parameters and are observed in the presence of Gaussian noise [6]. 
However, when the observables are nonlinear functions of the parameters, extended Kalman filters 
(EKF) are employed which basically linearize the functions locally around the current estimate. 
As a result, EKFs are not the optimum estimators in the absolute sense. 

For the current application, both third- and fourth-order EKFs were analyzed and simulated. Unlike 
the MLE, the EKF provides frequency estimates every 2 ms rather than an averaged estimate. 
Additionally, both phase and frequency rate estimates are provided. The transient response of the 
EKF was improved by including an additional exponential coefficient to weight past data. For the 
given trajectory, weighting the past 40 samples was found to provide the best performance. The 
details of the equations are all included in [5] along with the performance curves for different filters 
and parameters. Frequency threshold attained by the fourth-order EKF is about 24 dB-Hz with a 
3-Hz rms error. Although phase is of secondary importance in the current application, its threshold 
was only 0.5-dB higher at 24.5 dB-Hz with a 0.4-rad rms error. 


4.3 Cross-Product AFC (CPAFC) Loop 


A simple sub-optimal method of estimating the frequency of a sinusoidal signal embedded in noise 
and subject to severe dynamics is the CPAFC loop. Basically, frequency discrimination is obtained 
by employing a cross-product between the in-phase and quadrature signals, thus removing the 
phase. The error signal is then filtered and fed to a numerically controlled oscillator (NCO). Even 
though the CPAFC loop permits operation at low CNRs where a PLL is inoperative, it also increases 
the rms frequency error. 
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In the presence of time-varying Doppler, the loop bandwidth needs to be adjusted such that a 
balance is reached between noise immunity (obtained with small loop bandwidth) and dynamic 
tracking (typically requiring larger loop bandwidth). Both second- and third-order loop filters 
were tested for the given jerk [5]. A threshold of 24.7 dB-Hz was achieved with an 8-Hz loop 
bandwidth and second-order filter. However, the frequency estimation error was 31 Hz rms at 26- 
dB-Hz CNR. Using a third-order filter with loop bandwidth of 5.5 Hz can achieve a slightly lower 
threshold (approximately 0.3-dB lower) with almost no impact on rms frequency error. 


4.4 Digital Phase-Locked Loop (DPLL) 


Frequency estimation can be achieved with traditional DPLLs by processing the phase estimates. 
However, in regions where the DPLL is inoperative due to cycle slipping, no frequency information 
can be derived; hence, the frequency operation threshold is dominated by the phase behavior, which 
typically requires higher CNRs. 

A type-III DPLL was simulated [5] and it was found that a nominal loop bandwidth of 43 Hz 
minimized the probability of loss-of-lock for our trajectory. A threshold around 26 dB-Hz was 
achieved with an rms frequency error on the order of 26 Hz. 


4.5 Overlapping DFT AFC (ODAFC) Loop 


Frequency discrimination can be achieved in several ways, one of which is to perform a simple 
cross-product between the in-phase and the quadrature signals as we have seen before. However, 
it is straightforward to show [7] that the cross-product is equivalent to a 2-point overlapping DFT. 

A generalization would be to employ an N-point DFT where N can be optimized to offer good noise 
reduction while tracking dynamics. In order to maintain the loop update rate at 2 ms, overlapping 
DFTs are used, offering better performance. Analysis and performance of this improved loop 
are presented in [7]. For the 100-g/s trajectory, a 4-point FFT implementation with 10-Hz loop 
achieves a 22.5-dB-Hz CNR threshold and rms frequency error of 17 Hz. Using an 8-point DFT 
and increasing the loop bandwidth to 14 Hz reduces rms frequency error to 14 Hz with an increase 
in threshold CNR to 23 dB-Hz. 


4.6 Frequency Extended Kalman Filter (FEKF) 


Tracking algorithms that ignore phase information appear to have a lower threshold CNR (as 
frequency trackers) than their counterparts that also estimate phase. This is certainly the case 
when comparing the DPLL and the CPAFC loop considered above. However, improvement in 
threshold CNR is achieved at the expense of degraded performance in terms of frequency rms 
error. 
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This situation provides the basis for the frequency EKF which performs a cross-product on the data 
in order to remove the phase. This modified set of data is then used in an EKF with a reduced 
order. This system can also be thought of as a CPAFC loop with an optimized loop filter. The 
analysis of such an algorithm is presented in [8] along with computer simulations. A second-order 
FEKF was able to achieve a 22.5-dB-Hz CNR threshold for the given trajectory with a 41-Hz rms 
frequency error. This certainly outperforms the CPAFC loop described earlier in terms of lowest 
operating CNR. 


4.7 FFT-CPAFC Loop 


A block diagram of the CPAFC loop with FFT measurement correction as implemented by IEC 
is shown in Figure 4.7.1. This frequency estimator relies on a robust CPAFC loop to maintain 
lock, even in the presence of severe dynamics and additive noise. Although AFC loops can often 
maintain lock in extremely noisy environments, the rms estimation errors produced by these loops 
during dynamics are generally quite high. The central idea behind the IEC implementation is that 
if the instantaneous estimation errors could be estimated accurately and subsequently subtracted 
from the initial CPAFC estimates, then perhaps the total rms estimation error could be reduced. 
To this end, an FFT-based frequency estimator is employed to measure the residual error present 
in the in-phase and quadrature samples following the data-wipe operation, and prior to taking the 
cross-product. At this point, the 1-ms samples are converted to 2-ms samples within the CPAFC 
loop, and to 4-ms samples in the FFT estimator. Thus, the maximum frequency error possible 
within the CPAFC loop without loss-of-lock is ±250 Hz, while the maximum range of the FFT 
estimator is ±125 Hz. The output of the CPAFC drives the NCO, whose output closes the loop. 
An averaged version of the CPAFC output is also available every 100 ms for combining with the 
FFT estimates. 

The FFT estimator operates as follows: every 100 ms, 25 4-ms samples are collected, 7 zeros 
appended, and a 32-point FFT performed. A sinusoidal component in the samples gives rise to an 
impulse in the FFT output near the proper frequency, with a resolution of roughly 8 Hz (250/32). 
Thus, the magnitude of the FFT is computed, and the frequency corresponding to the location 
of the peak declared to be a rough estimate of the AFC frequency error. A refined estimate is 
obtained by performing 3-point quadratic interpolation. This refined estimate is then combined 
with the averaged CPAFC loop estimate synchronously every 100 ms to obtain the final frequency 
estimate. 

The perfomance of this hybrid frequency estimator has been evaluated by means of simulation, 
nonlinear analysis, and hardware experiments. The results of the hardware evaluation are presented 
in Section 6.1.4, while the simulation resuts and the analysis are presented in Appendix A. 
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Figure 4.7.1: FFT-CPAFC block diagram (figure courtesy IEC) 
















4.8 Adaptive Least Squares (ALS) 


Several algorithms were investigated under the collective title of Adaptive Least Squares. These 
algorithms are recursive in nature and require higher computational complexity than the previously 
described methods, but have potentially better performance. An ideal application for these algo- 
rithms is the post-mission processing of data. The following paragraphs describe these algorithms 
briefly, with more detailed descriptions in [9] and [10]. 


4.8.1 Fast Frequency Acquisition 

This algorithm provides for fast carrier acquisition in a data-wipe environment. The received signal 
is first demodulated by a carrier reference signal of known frequency and phase and by its 90-degree 
phase shifted version. The basic equations for the measurement express the noisy samples as a 
truncated (nth order) series involving the unknown frequency and phase, the sampling times, and 
the noise samples. The equations [9] are put in linear form involving a matrix of parameters that 
depends on the frequency and phase, and from which the frequency and phase can be determined. 

The algorithm obtains the matrix of parameters from a sequence of N pairs of measurements using 
least squares estimation. This requires inversion of an N x N matrix that contains terms dependent 
on the sampling times and that can be precomputed and then multiplied by a matrix formed 
from the measurements. The matrix to be inverted has a specific structure that makes it possible 
to use a rapid algorithm for the solution. As a result, the desired matrix of parameters can be 
obtained in about 6N\og 2 N operations, much fewer than the number of operations required for the 
“brute-force” calculations of a general matrix equation of similar form. 


4.8.2 Differential Sampling for Fast Acquisition 

This is an extension of the previous algorithm using a differential signal model along with appropri- 
ate sampling techniques. The algorithm is recursive in measurements and thus the computational 
requirements increase only linearly with the number of measurements. 

The dimension of the state vector in the proposed algorithm does not depend on the number 
of measurements and is quite small, typically around four. This is an advantage compared to 
the previous algorithm where the dimension of the state vector increases monotonically with the 
product of frequency uncertainty and observation period. Further details are provided in [10]. 


4.9 Signal Acquisition Algorithms 

The initial acquisition of a GPS signal consists of two phases: first the individual tracking channels 
estimate pseudorange and range rate to the respective satellites, then, when sufficient channels 
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have acquired signals, a state (position, velocity, and time) solution is derived. Results presented 
here are concentrated on the first phase, channel acquisition. 

The problem can be viewed as a two-dimensional search, as shown in Figure 4.9.1. The unknown 
pseudorange is plotted on the horizontal axis, while the unknown range rate is plotted on the vertical 
axis. Pseudorange is obtained by searching for the appropriate C/A code lag that maximizes the 
code correlator output. Since the C/A code has 1023 lags, this search requires 1023 correlations 
with one-chip spacing or 2046 correlations with half-chip spacing. (The latter allows operation at 
lower CNR.) The search can be performed serially, i.e., sequential dwells for each possible code 
lag, using a single correlator, or in parallel, i.e., obtain all correlations during a single dwell by use 
of a bank of correlators. Parallel searching is feasible with today’s technology using VLSI. As an 
example, the TPS uses two VLSI chips to perform all 1023 correlations. 

The second search is for range rate. Unknown range rate is due to two factors: uncertainty as to the 
center frequency, e.g., crystal frequency offset, and unknown Doppler due to the mutual dynamics 
of the vehicle and GPS satellite. Here also the search can be performed serially or in parallel. 
In the sequential search, used in the TPS, each step covers a frequency range that is inversely 
proportional to the integration time. The TPS integration time is 1 ms, hence the frequency 
coverage of a single search step is approximately from -500 Hz to 500 Hz, with significant CNR 
degradation at higher frequencies. In each step of the serial search, the channel’s carrier NCO is 
positioned to a different center frequency so the search is for a range around that center frequency. 
In a parallel search, correlator outputs are sampled at a much higher rate and processed via FFT 
processors, hence a significantly higher frequency range can be covered. If total integration time 
remains fixed, then sampling the correlator output at 32 times per integration time will result in 
a searched frequency range that is 32 times larger than that covered by a serial frequency search. 
However, this improvement is at the cost of additional hardware complexity. 

Algorithms for fully parallel searches are presented and analyzed in [11] and [12]. A summary of 
results is presented in Figure 4.9.2. The figure illustrates channel acquisition time as function of 
Doppler uncertainty, for various CNR. We observe that for RAJPO’s specification of a Doppler 
uncertainty of 6150 ft/s and a CNR of 28 dB-Hz (using the 3-dB loss assumption), the channel 
acquisition time is less than 0.5 s. Even when adding the settling time of the navigation filter, this 
approach can save approximately 1 s from RAJPO’s TTFF specification. 

The comparison of rms frequency error is somewhat unfair to the phase tracking algorithms. When 
the estimates are averaged over long time intervals, the rms frequency error becomes much better 
for the phase tracking algorithms. This is because the improvement is linear in time when phase 
tracking, but is proportional to the square root of time when not phase tracking. 


4.10 Performance Comparison via Simulations 


We conclude this chapter by comparing the different frequency tracking algorithms for the “data 
wipe” environment, using the trajectory described earlier in Section 3.3. 

Table 4.1 summarizes the simulation results for the frequency tracking schemes. The comparison 
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Figure 4.9.2: Acquisition time versus Doppler uncertainty 
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Table 4.1: Performance Comparison via Simulations 


Algorithm 

Threshold 

(dB-Hz) 

RMS Er 
23 dB-Hz 

ror (Hz) 

26 dB-Hz 

Phase 

Estimate 

Freq. Rate 
Estimate 

Complexity 

ODAFC 

22.5 

17.0 

9.0 

No 

No 

3 

FEKF 

22.5 

36.0 

22.5 

No 

Yes 

4 

ML 

23.0 

7.0 

1.0 

No 

Yes 

7 

EKF 

23.9 

3.5 

2.2 

Yes 

Yes 

6 

FFT-CPAFC 

24.0 

38.0 

11.0 

No 

No 

5 

CPAFC 

24.7 

60.0 

31.0 

No 

Yes 

2 

DPLL 

25.7 

20.0 

12.0 

Yes 

No 

1 


includes operating threshold (defined as the CNR at which the probability of frequency loss-of-lock 
is 10 percent) and rms frequency error at a low CNR (23 dB-Hz) and at a relatively higher CNR 
(26 dB-Hz). Other valuable information, such as whether phase and/or frequency are estimated 
and algorithm complexity (1 denoting the least demanding algorithm), are listed. The algorithms 
are listed in ascending order by threshold CNR. Lowest threshold CNR, 22.5 dB-Hz, is achieved 
by the ODAFC loop and the FEKF. MLE and EKF algorithms have the smallest rms frequency 
errors (a few Hertz) with the MLE exibiting almost 1-dB lower threshold CNR than the EKF at 
23.0 dB-Hz. Depending on whether the operating threshold or the rms frequency error is of prime 
importance in a specific application, a suitable algorithm can be selected from Table 4.1 to satisfy 
performance requirements and available computational power. 
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Chapter 5 


Description of Experimental Setup 


The setup of the GPS high-dynamic receiver is considered and presented in this chapter. The 
breadboard Advanced Receiver (ARX) is described first since it constitutes the testbed for the 
GPS demonstration. Software as well as hardware modifications to the ARX are then discussed 
with emphasis on signal generation, as the latter contains the dynamics of interest. 

The ARX is used as a testbed because it represents hardware typical of modern GPS receivers, 
specifically the TPS. The ARX obtains digital I and Q samples at a high rate (5 MHz for the ARX 
versus 2 MHz for the TPS), integrates the samples to reduce the data rate (1-ms integration for 
the TPS versus variable integration time in the ARX), and performs most of the remaining pro- 
cessing in software (Intel 80286 in the ARX versus Motorola 68000 in the TPS). Hence, algorithms 
demonstrated on the ARX have practical application to the TPS. 


5.1 Advanced Receiver Block Diagram 

The ARX is a hybrid analog/digital receiver which has been designed for use in NASA’s Deep Space 
Network (DSN). It uses intermediate frequency sampling and DPLL’s to perform carrier tracking, 
subcarrier tracking, and symbol synchronization [13]. 

Figure 5.1.1 shows the ARX signal processing diagram. The receiver uses the open loop IF signal 
at 53 MHz from the RF front end. This signal is then passed through a total power automatic gain 
control (AGC) circuit before entering the carrier loop. The reference signal is produced by mixing 
a constant 46.25-MHz signal with the output of the carrier NCO (nominal 1.875 MHz) frequency. 
The output from this mixing passes through a lowpass filter and is digitized. The clock driving the 
A/D is derived from the symbol synchronization loop NCO and has a nominal value of 19.75 MHz. 
This technique ensures a fixed number of samples per symbol when the signal is BPSK modulated. 
The digitized 5-MHz signal (sampled at roughly 20 MHz) is then split into in-phase and quadrature 
components (each with a 10-MHz rate). Conversion to baseband is accomplished digitally using 
a 5-MHz signal phase locked to the sampling clock. The 10-MHz rates are then lowered to the 
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Figure 5.1.1: ARX signal processing diagram 

























appropriate loop update rates using digital accumulators to perform carrier tracking, subcarrier 
tracking, and symbol synchronization. The various loop filters can also be set independently in 
software with different parameters, including loop order and loop bandwidth. 


5.2 Modifications to the ARX 


In the current application, the received signal consists of a pure sine wave with a time-varying 
frequency driven by the prescribed trajectory. The amplitude of the tone is set such that the 
required CNR is achieved with fixed noise power. Subcarrier and symbol synchronization are not 
used since they are not required for tracking the frequency of an unmodulated tone. 

The required software modifications constituted a significant part of the work, which included 
algorithm implementation, system programming, and performance evaluation. Hardware changes 
were required mainly in the signal generation subsystem for the purpose of generating the dynamics. 


5.2.1 Hardware Modifications 

Modifications in the system’s hardware were essential due to the nature of the experiment. The 
A/D clock was set to a fixed 19.5 MHz derived directly from a frequency synthesizer. As for the 
signal generation, the software-generated trajectory was driving an NCO with a 1.5-MHz nominal 
frequency, its output was mixed with a fixed 18-MHz reference, band-pass filtered, and then up- 
converted again to 53 MHz using a fixed 33.52-MHz reference. The center frequency of 53.02 MHz 
was chosen to avoid a bias observed at 53 MHz. This bias will later be described and characterized 
in more detail. The new configuration is depicted in Figure 5. 2. 1.1 with the signal generation path 
shown in Figure 5. 2. 1.2. 


5.2.2 Software Modifications 

The Advanced Receiver software program (ARX) is a menu-driven system that allows the operator 
to set parameters such as bandwidth, symbol-to-noise ratio (SNR), carrier and subcarrier update 
rates, etc. It also allows the operator to choose from several possible displays for the CRT monitor 
and, upon operator request, will log statistical data on a disk. The ARX, upon request, performs 
either preacquisition or tracking. 

The ARX is divided into sixteen sections called modules. Eight of these sixteen modules had 
to be modified to accommodate the differences in the requirements of the ARX and the Global 
Positioning Satellite (GPS) demonstration. The eight modules that needed modification are: 


1. ADRCF — the foreground module that acquires the in-phase and quadrature (I and Q) samples 
from the A/D convertors, solves the loop filter equations, drives the NCO (see Figure 5. 2. 1.1) 
and acquires data for displaying or logging. 
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2. ARPRINT — the background module that performs the data reduction, writes to the CRT 
monitor and logs to the disk. 

3. ARFAST — the module that performs all of the communication with the array processor board 
(AP-4) for FFT processing as part of frequency error estimation. 

4. ARPCALC — the module that calculates all of the parameters necessary to initialize the hard- 
ware (i.e., the count registers that control I and Q update rates), the loop filter coefficients, 
and the gains and constants needed for data reduction. 

5. ARINITL — the module that sets the default values to parameters and clears memory areas 
where required. 

6. ARPENTR — the module that accepts and translates the operator inputs for setting param- 
eters and controlling the execution of the software. 

7. ARMENUS — the module with the monitor displays of the menus. 

8. ARGLOB — the insert module that defines all of the globally used parameters. 

The changes to each of the eight modules of the ARX to effectively simulate the GPS are described 
in Appendix C. 
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Figure 5.2. 1.1: ARX hardware block diagram 
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Figure 5.2. 1.2: Trajectory generation hardware 


27 







Chapter 6 


Algorithm Evaluation 


The demonstration tests are divided into two categories: static tests, which are conducted with a 
constant signal frequency at a specific CNR, and dynamic tests, which also include a simulated tra- 
jectory. Not all algorithms described in Chapter 4 have been demonstrated. The DPLL, CPAFC 
loop, EKF, and FFT-CPAFC loop algorithms have all been tested with and without dynamics. 
The MLE algorithm was not tested due to schedule limitations. The adaptive least squares algo- 
rithms could not be tested as they require different setup conditions such as data modulation. The 
remaining two algorithms, the ODAFC loop and the FEKF, were not part of the demonstration 
plan. 

In the static tests, only rms frequency error is measured for different loop or algorithm parameters. 
For dynamic tests, results also include estimates of the probability of losing frequency lock. The 
latter is defined to have occurred when the instantaneous frequency error exceeds in magnitude 
half the Nyquist sampling rate. In this application, the loop update rate is set to 500 Hz and 
loss-of-lock is declared when the instantaneous frequency error exceeds ±250 Hz. 

Figure 6.1 depicts the signal processing of the ARX relevant to the GPS demonstration. The 
input signal is mixed to an IF frequency of 5 MHz, then is bandpass filtered and sampled. All 
subsequent signal processing is digital, including removing the 5-MHz IF and generating in-phase 
(I) and quadrature (Q) baseband samples, which are integrated over 2 ms. The loop is closed via a 
high resolution (32-bit) NCO that is clocked at 10 MHz. Each tested algorithm can use either both 
I and Q samples or only the Q samples (as with the DPLL), process these samples, and output a 
frequency estimate to the NCO every 2 ms. 

The input to each algorithm consists of a sequence of I and Q samples normalized to unity ampli- 
tude. This can be accomplished because the amplitude of the received signal is assumed known to 
the receiver. An equivalent model for the signal path to determine the various loop gains is shown 
in Figure 6.2. Let A denote the amplitude of the received signal. The reference signal has unit am- 
plitude which renders the analog mixer’s output to an amplitude equal to A/2. The digital complex 
mixers following the demultiplexer do not alter the amplitude. After the mixer is a two-sample 
accumulator followed by an Mi -sample digital summer to reduce the processing rate. In this case, 
M\ is equal to 19500 and it reduces the rate from 9.75 MHz (half of the 19.5-MHz sampling rate) 
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10 MHz 


Figure 6.1: Demonstration data flow 


to the update rate of 500 Hz. The accumulator output is 35 bits long. Since the Intel 80286 CPU 
used in the system uses 32-bit numbers, a division by eight shown in Figure 6.2 is used for proper 
scaling. In order to normalize the samples, a final division by ( AM\/8 ) is performed that renders 
the value of the amplitude of I and Q to unity. All algorithms process these samples to derive the 
frequency estimate. The NCO accepts a 32-bit integer and hence the frequency estimate needs to 
be multiplied by K /v = 2 32 /10 7 to convert to the NCO units. 

Since the hardware differs slightly from the simulation model, one expects the hardware and simu- 
lation tests to be slightly different. The following possible discrepancy sources have been identified: 

1. The actual amplitude of the 2-ms sample is not constant, as assumed in the simulation, but is 
proportional to sin( Au/T 3 Mi/2)/Misin( Au>T,/2) where A lj is the frequency error at the input 
to the code correlator and T s is the sampling period. This primarily impacts the probability 
of loss-of-lock of the algorithms. 

2. The calibration of the hardware signal generator, evaluated in Appendix B, is on the order 
of 0.5 dB. This is a CNR-dependent calibration. 

3. The nature of the noise in the actual system slightly deviates from the AWGN assumption. 
For example, there is a slight asymmetry in the passband filter. 

4. The timing of loading the NCOs in the hardware system causes an extra transport lag of 
0.225 ms compared to the software simulation. This effectively changes the loop bandwidth 
and causes a more sluggish response to dynamics. 


These differences and others contribute to the difference between hardware and software results, 
but are estimated to have a total effect of less than 2 dB in threshold CNR. 
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Figure 6.2: Equivalent loop gain model 


6.1 Static Test Results 


The static tests do not require any hardware modification since the frequency of the generated 
signal is constant. Only a 19.5-MHz clock is needed to drive the sampler and that was supplied 
from a frequency synthesizer. Implementation of the algorithms is accomplished in software except 
for those using FFT, which were computed using a dedicated array processor board. 


6.1.1 DPLL 

The DPLL was already implemented in the ARX and therefore the tests were straightforward to 
run. The algorithm consists of a third-order filter which uses only the Q samples and is given by 

*W-o> + l3=r + jr^j! < M > 


where 




G 3 = 


krd 3 

~ 
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Figure 6. 1.1.1: DPLL static test results 


^ _ 4 BcT(r — k) 
r(r - k + 1) 

Here, T denotes the 2-ms update time and r and k are filter parameters. Performance of the type- 
III loop is shown in Figure 6. 1.1.1 for different loop bandwidths and k = 0.5, r = 3. Bandwidths 
tested correspond to 20, 40, and 60 Hz in anticipation of the optimum 43-Hz loop bandwidth to be 
used in conjunction with the trajectory. 

As shown in the figure, the simulation and measured points are in agreement to within 0.5 dB. 
Each measured point was obtained using 140,000 samples, while the simulation points were carried 
over 50,000 samples. 
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6.1.2 CPAFC 


The CPAFC frequency discriminator uses both the I and Q samples in tracking the received fre- 
quency. The discriminator output is represented by 


V(k) = I(k - l)Q(k) - Q(k - l)I(k) 


( 6 . 2 ) 


where k denotes discrete time. This is followed by a filter with transfer function 


F{z) 


G i G-j 

1 - 2- 1 (1 - Z - 1 ) 2 


(6.3) 


where G \ and G 2 were given previously. Test results are shown in Figure 6. 1.2.1 for a second-order 
loop (r = 2) with an update rate of 500 Hz. The bandwidths tested correspond to 5, 7, and 10 Hz in 
anticipation of the optimum 7.5-Hz loop bandwidth to be used in conjunction with the trajectory. 

As shown in Figure 6. 1.2.1, simulation and measured points are in agreement to within 0.5 dB. 
Each measured point was obtained using 200,000 samples while the simulation points were carried 
over 500,000 samples. 

The tests were carried out at 53.02 MHz, a 20-kHz deviation from the nominal center. The deviation 
was necessary to reduce the effects of a bias, nonzero mean frequency error, encountered at 53 MHz. 
The cause of the bias has not been determined, but the bias was both CNR and center frequency 
dependent. Figure 6. 1.2.2 depicts the bias at CNR = 30 dB-Hz for a 5-Hz loop bandwidth for both 
53- and 53.02-MHz center frequencies. As shown, the estimation error is biased at 53 MHz, and 
for all practical purposes, unbiased at 53.02 MHz. 


6.1.3 EKF 

The EKF was implemented with steady-state gains to reduce the computational burden of com- 
puting and updating several vectors and matrices. For a fourth-order filter, the algorithm is given 

by 


x\(k) = g(k - 1) + kin(k) (6.4) 

T 2 

x 2 (k ) = x 2 (k - 1) + Tx 3 (k - 1) + — x 4 (k - 1) + k 2 fi(k) (6.5) 

i 3 (A:) = xs(k - 1) + Tx^(k - 1) + kzn(k) (6.6) 
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Figure 6. 1.2.1: CPAFC loop static test results 
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x 4 (k) = x 4 (k — 1) + k 4 fi(k) 


(6.7) 


g(k) = xi(fc) + Tx 2 (fc) + 


T 2 

2 


2^3 

*3(fc) + -g-au(fc) 


( 6 . 8 ) 


where /x(fc) is the quadrature sample normalized by ^±y/2P c . The frequency fed to the NCO is 
derived from g(k) as 


f(k) = 


</(fc) - g(k - 1) 
27T 


(6.9) 


where k\, fc 2 , £3, and k 4 constitute the first column of the steady-state error covariance matrix 
Y given in [8]. Setting k 4 to zero results in a third-order filter. It is worthwhile noting that the 
implementation corresponds to an approximation of the steady-state EKF in the sense that the 
signal controlling the NCO was set to be the difference between two consecutive phase estimates, 
whereas strictly speaking, the EKF requires mixing the received I and Q samples with an exact 
phase. The loss due to that approximation is shown in Figure 6. 1.3.1 and is about 1.7 dB. For 
that reason, the measured threshold in the presence of dynamics will be at least 1.7 dB worse than 
previously predicted in [8]. 

In Figure 6. 1.3.2, the simulation and measured points are compared for both filter orders and are 
found to be in agreement. The phase, frequency, frequency rate, and the derivative of the frequency 
rate are given by xi(fc), x 2 (k), X3 (fc), and x 4 (k), respectively in units of rad, rad/s, etc. 


6.1.4 FFT-CPAFC 

The FFT-CPAFC loop was implemented in the ARX using an AP-4 digital signal processing board 
manufactured by DSP Corporation. The implementation followed the description given in Section 
4.7, except that the minimal sampling rate was 2 ms throughout, instead of the 1-ms primary 
sampling rate employed by IEC. Two consecutive I and Q samples were summed prior to FFT 
processing, yielding an FFT frequency range of ±125 Hz. The 25 in-phase samples were pre- 
multiplied by an alternating sequence of plus and minus signs in order to shift the zero frequency 
point to the center of the frequency range. The index of the greatest magnitude in the 32-point 
FFT output was declared to be the rough estimate of CPAFC loop frequency error. Three-point 
quadratic interpolation was performed on the rough estimate to obtain the final estimate, which 
was then subtracted from the averaged CPAFC frequency estimates. In this manner, a refined 
frequency estimate was obtained every 100 ms. 

Static tests used 100 independent runs at each CNR. RMS frequency error is shown in Figure 
6. 1.4.1, indicating that hardware results slightly outperform the simulations. 
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Figure 6. 1.3.1: Implementation loss in the EKF 
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Figure 6. 1.3.2: EKF static test results 
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Figure 6. 1.4.1: FFT-CPAFC loop static test results 




6.2 Dynamic Test Results 


Generation of the trajectory for the dynamic tests was accomplished by coding in software a subrou- 
tine that generates the frequency of the Doppler at 2-ms intervals, adds it to the 1.5-MHz nominal 
frequency of NC01 (Figure 5.2. 1.1), converts it to a 32-bit integer, and finally outputs it to the 
NCO. The test trajectory is described in Chapter 3. 


6.2.1 DPLL 

DPLL dynamic tests were conducted for different CNRs, using a 43-Hz loop bandwidth. Results 
are shown in the attached figures for a third-order loop with an update rate of 500 Hz ( Bl = 43 
Hz, r = 2, k = 0.5). 

Transient frequency error versus time is plotted at very high CNR (73 dB-Hz in Figure 6. 2. 1.1 and 
60 dB-Hz in Figure 6.2. 1.2), showing steady-state error due to jerk at about 0.8 Hz. Measurement 
data are in agreement with the simulations (carried out without any noise) shown in solid line. 
Measured and simulated frequency rms error and probability of loss-of-lock are shown in Figures 
6.2. 1.3 and 6. 2. 1.4 versus CNR. The measured threshold (CNR at which frequency looses lock with 
0.1 probability) occurred at 26 dB-Hz, while the simulated threshold occured at 25 dB-Hz. As for 
the rms values, they are in agreement to within 0.8 dB, as shown in Figure 6.2. 1.3. 


6.2.2 CPAFC 

Tests were run with a second-order filter using an 8-Hz loop bandwidth and r = 2. For the given 
trajectory and filter, these parameters were found by simulation to minimize the CNR at which 
loss-of-lock occurred with a 10 percent probability. 

Figure 6.2. 2.1 presents transient frequency error versus time at very high CNR (60 dB-Hz) to show 
the steady-state error due to jerk, which is about 14.5 Hz. Measurement and simulation data are 
in close agreement. Frequency rms error and the probability of loss-of-lock versus CNR are shown 
in Figures 6.2. 2.2 and 6. 2.2. 3, respectively. Measured rms error is in close agreement with the 
simulations but the measured threshold is about 1.4 dB higher at roughly 24 dB-Hz. 


6.2.3 EKF 

Both the third- and fourth-order EKFs were tested in the presence of dynamics (the fourth-order 
EKF has been shown by simulation to possess a lower CNR threshold). 

Figure 6.2. 3.1 shows frequency error versus time at 60 dB-Hz, with maximum error due to jerk of 
4.2 Hz for a third-order EKF. The measured points agree well with the simulations (also carried 
out at 60 dB-Hz) shown in dotted line. Frequency rms error and the probability of loss-of-lock 
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Figure 6. 2. 1.3: DPLL dynamic test results (RMS frequency errors versus CNR) 







TIME, s 


Figure 6.2.2.1: CPAFC loop transient response (CNR = 60 dB-Hz) 
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RMS FREQUENCY ERROR, Hz 



Figure 6. 2. 2.2: CPAFC loop dynamic test results (RMS frequency error versus CNR) 




PROBABILITY OF LOSS-OF LOCK 




versus CNR are shown in Figures 6.2.3.2 and 6. 2. 3.3. The measured rms error is in close agreement 
with the simulations but the measured threshold of the approximate fourth-order EKF is about 1 
dB higher at roughly 26 dB-Hz. 


6.2.4 FFT-CPAFC 

Dynamic tests were carried out using the IEC trajectory, described in Chapter 3. Sample statistics 
were obtained at various CNRs between 20 and 30 dB-Hz. Results are shown in Figure 6.2.4. 1, 
where good agreement between experiment and theory is apparent. 

The probability of loss-of-lock was also estimated for the dynamic case. Loss-of-lock occurs when 
the instantaneous CPAFC error exceeds ±250 Hz. Estimated loss-of-lock probabilities are shown in 
Figure 6.2.4. 2, along with the simulated results. The simulation indicates that loss-of-lock becomes 
a problem near 22 dB-Hz, while for the experiment, this threshold is closer to 24 dB-Hz. This is 
attributed to amplitude attenuation that is unmodeled in the simulations (see Section 6.3). 

Performance of the third-order CPAFC (and other third-order loops) is very sensitive to loop 
bandwidth. Since maximum transient error due to jerk is roughly inversely proportional to the 
third power of the bandwidth, small changes in bandwidth result in large changes in peak error. 
Since the CPAFC loses lock when the instantaneous frequency error is over 250 Hz, 125 Hz for the 
FFT, small changes in bandwidth or input dynamics can result in major changes in the performance 
curves. An example is shown in Figure 6. 2.4.3, which is a plot of the instantaneous frequency error 
of the CPAFC (solid curve) when 150-g/s jerk is applied, along with the FFT estimate (long 
dashes), and the combined frequency error (short dashes). At roughly 1.5 seconds into the test, 
the CPAFC error exceeds 125 Hz, which is the upper limit of the FFT estimator. This results in 
aliasing, whereby the FFT peak suddenly shifts to -125 Hz, yielding an error of roughly 250 Hz in 
the final estimate. This phenomenon can lead to catastrophic performance degradation, as shown 
in Figure 6. 2.4.4. The situation is easily remedied by increasing the loop bandwidth somewhat (in 
this case to 7 Hz) in order to reduce the CPAFC’s dynamic excursions. 


6.3 Comparison and Discussion 


Performance of the tested algorithms are compared in this section for static and dynamic tests. 
Performance measures consist of rms frequency error and probability of loss-of-frequency-lock in 
the presence of dynamics only. Comparisons are made on the basis of measured data, without the 
application of the calibration curve of Appendix B. 

Table 6.1 shows results of rms frequency error comparisons for dynamic tests. Hardware demon- 
stration and simulation results for the various algorithms are in close agreement, within 5 percent, 
as shown in this table and in previous figures. The fourth-order EKF performed with the lowest 
rms frequency error, 3.2 Hz at 26 dB-Hz, but its operating threshold was at 25.8 dB-Hz. 

Results for threshold CNR, shown in Table 6.2, exhibit an average 1.4-dB difference between 
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Figure 6. 2. 3.2: EKF dynamic test results (RMS frequency error versus CNR) 
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Figure 6.2.4.1: FFT-CPAFC dynamic test results (RMS frequency error versus CNR) 




CNR.dB-Hz 


Figure 6. 2.4.2: FFT-CPAFC dynamic test results (probability of loss-of-lock versus CNR) 
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TIME, s 

Figure 6.2.4.3: FFT-CPAFC transient response (CNR = 30 dB-Hz) 





hardware and simulation results, with hardware threshold exceeding simulated (software) values. 
(Adding calibration data raises the difference to 1.9 dB.) The lowest measured threshold CNE, 
23.8 dB-Hz, was achieved by the CPAFC loop, while the highest was in the third-order EKF. The 
DPLL CNE threshold of 26 dB-Hz is 2.2 dB higher than the CPAFC loop. 

The fact that hardware and software agree well in rms frequency error while showing inferior 
hardware performance in CNE threshold can be traced to the effect of the amplitude attenuation 
of input samples due to the digital accumulators that was not modeled in the software simulations. 
This effect is significant only when the instantaneous frequency error is a large fraction of 250 Hz. 
To validate this hypothesis we re-tested a third-order CPAFC with a bandwidth of 5.2 Hz, using 
the 70-g/s trajectory. Figure 6.3.1. a shows the rms frequency error when the amplitude attenuation 
effect is present and absent, with no significant difference in rms frequency error. Figure 6.3. l.b 
shows that, for the same configuration, amplitude attenuation degrades CNE threshold by 1-2 dB. 
This attenuation affects the various algorithms in different ways, depending on the frequency error 
during the maneuvering. For instance, the DPLL has a 0.8-dB difference between simulated and 
measured CNE threshold because of its smaller frequency error in the presence of jerk, as shown in 
Figure 6. 2. 1.1. Depending on the algorithm in question, it is believed that the effect of amplitude 
attenuation is 0.8 dB to 1.8 dB in CNE threshold. 

Simulated threshold CNEs in Table 6.2 are also slightly different than those of Table 4.1 because 
of the different trajectories involved. The initial trajectory of Table 4.1 included Doppler rate at 
the start of the simulations, while the trajectory of Table 6.1 contains a constant initial frequency. 

In conclusion, the tested algorithms performed within 0.5 dB of simulations in rms frequency error 
and within 1.5 dB of simulations in loss-of-lock threshold. The CPAFC loop exibited the lowest 
threshold and the EKF achieved the smallest rms frequency error. 
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CPAFC LOOP (THIRD-ORDER) 

(B L = 5.2 Hz, r = 2, k = 0.5) 

WITHOUT AMPLITUDE ATTENUATION 

WITH AMPLITUDE ATTENUATION 

70 g/s TRAJECTORY 



CPAFC LOOP (THIRD-ORDER) 

(B L = 5.2 Hz, r = 2, k = 0.5) 

WITHOUT AMPLITUDE ATTENUATION 

WITH AMPLITUDE ATTENUATION 

70 g/s TRAJECTORY 


CNR, dB-Hz 


Figure 6.3.1: CPAFC - impact of amplitude attenuation on (a) RMS frequency error; and 
b) loss-of-lock probability 




Table 6.1: Dynamic Test Results - RMS Frequency Error (Hz) 


Algorithm 

Trajectory 

(25 dB-Hz) 

(30 dB-Hz) 



Simulated 

Measured 

Simulated 


DPLL 


17.5 

18.0 

9.6 

9.8 

CPAFC 


19.9 

21.6 

9.7 


EKF (3) * 

1 

5.9 

6.1 

3.7 

3.7 

EKF (4) * 


3.5 

4.2 

2.4 

2.4 

FFT-CPAFC 

mmSSBi 

10.0 

10.4 

0.98 

0.93 

* Tested at 29 dB-Hz rather than 30 dB-Hz 


Table 6.2: Dynamic Test Results - Threshold CNR 


Algorithm 

Trajectory 

Threshold CNR (dB-Hz) 

Difference 



Simulated 

Measured 

(dB-Hz) 

DPLL 


25.2 

26.0 

0.8 

CPAFC 


22.6 

23.8 

1.2 

EKF (3) 


26.0 

27.6 

1.6 

EKF (4) 


24.8 

25.8 

1.0 

FFT-CPAFC 


21.5 

23.6 

2.1 
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Chapter 7 


Conclusions 


Frequency tracking at low CNRs and in the presence of high dynamics has been the goal of this 
research effort. Several novel schemes have been introduced, analyzed, and tested in order to 
validate their performances. These algorithms meet RAJPO’s requirements for operation of the 
TPS, and offer potential improvements to the existing techniques. 
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Appendix A 


Description and Analysis of the 
FFT-CPAFC Loop 


This appendix is devoted to the description, evaluation, and nonlinear analysis of the FFT-CPAFC 
loop used in the TPS. Simulation results are presented in Section A.l, and the nonlinear perfor- 
mance analysis is the subject of Section A. 2. 


A.l Simulation Results 


The simulation follows the block diagram of Figure 4.7.1, except that 2-ms update rates are used in 
the CPAFC loop instead of IEC’s multi-rate implementation, and data- wipe is assumed. A type-III 
CPAFC loop was simulated with a single-sided loop bandwidth of 5.2 Hz, and with loop parameters 
r = 2, k — 0.5. Estimator performance is evaluated using the 100-g/s trajectory described in Section 
3.3, and the effects of additive noise are incorporated by adding independent Gaussian noise samples 
to the signal samples in the I and Q channels. The noise variance is adjusted to yield the desired 
CNR. Instantaneous frequency estimation errors were recorded at the output of the CPAFC loop 
(point B in Figure 4.7.1), the output of the FFT error estimator (point C), and the combined 
estimator output (point D). A characteristic sample sequence of instantaneous estimation errors is 
shown in Figure A.l.l.a at a CNR of 30 dB-Hz. Note that at this CNR the FFT estimates the 
CPAFC loop error accurately, yielding a combined rms error of only a few Hertz rms. At lower 
CNRs, the FFT estimator begins to suffer from outliers, as can be seen in Figure A.l.l.b ( CNR 
= 25 dB-Hz), resulting in degraded performance. At still lower CNR, the CPAFC loop begins to 
lose lock with high probability. An example of loss-of-lock behavior is shown in Figure A.l.l.c. 
Loss-of-lock becomes a serious problem below about 22 dB-Hz for this trajectory. At very high 
CNRs, the dominant error component in the combined estimate is due to dynamics, as shown in 
Figure A.l.l.d. Therefore, we expect the rms estimation error for the entire trajectory to become 
independent of CNR for CNRs well above 30 dB-Hz. 

The combined frequency estimation error for the entire trajectory was estimated at CNRs ranging 
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Figure A.l.l.a: Transient response to 100-g/s trajectory (CNR = 30 dB-Hz) 
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Figure A.l.l.b: Transient response to 100-g/s trajectory (CNR = 25 dB-Hz) 
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Figure A.l.l.c: Transient response to 100-g/s trajectory (CNR = 22 dB-Hz) 
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Figure A.l.l.d: Transient response to 



from 20 to 50 dB-Hz, using sample statistics obtained from 100 independent simulations at each 
CNR. The results are shown in Figure A. 1.2, where we also include simulation results for a static 
trajectory (zero dynamics), and the Cramer-Rao lower bound (also zero dynamics) for comparison. 
The discrepancy between the lower bound and the actual performance of the IEC estimator can be 
attributed partly to the presence of dynamics, and partly to the sub-optimality of the estimator 
structure. Estimated loss-of-lock probabilities are shown in Figure A. 1.3, which indicates that this 
algorithm should be able to maintain lock above 22 dB-Hz most of the time when tracking the IEC 
trajectory. 

When these tests were repeated with the ARX hardware, significantly worse performance for the 
100-g/s trajectory was recorded. The difference is mostly due to the amplitude attenuation effect, 
described in Section 6, that was not included in the simulations. In other words, the 100-g/s 
dynamics are in a zone where performance is highly sensitive to small changes in loop model or 
parameters. Repeated tests at 70 g/s show that hardware and software agree well. RAJPO’s 
specified dynamics, 50 g/s for 1 s, are well within the safe area. 


A. 2 Nonlinear Analysis 


This section analyzes the FFT-CPAFC loop in detail, specifically near threshold CNR. The analysis 
validates simulation results. 1 Referring to Figure A. 3, rms frequency errors at points A, B, C, and 
D are computed for a tone in the presence of AWGN, with and without dynamics. 


A. 2.1 Static Case 

Using the results of [7], if the loop noise bandwidth Bl < 10 Hz and CNR < 50 dB-Hz (but high 
enough to justify a linear analysis), the variance of the linearized loop is given by: 



Hz 2 


(A.l) 


where T a is the loop filter update time interval. Several techniques are available to predict the 
performance at lower CNRs [14]. The quasi-linearization method is especially attractive here, 
because of its simplicity. Omitting the details, since they are similar to those of [15], the quasi- 
linear performance is given by: 



(A.2) 


a = 


CNR 2 T s 


Bl 


(A.3) 


1 V. Vilnrotter and S. Aguirre, “Evaluation of IEC’s Frequency Estimation Algorithm,” Jet Propulsion Laboratory, 

Interoffice Memorandum (internal document), July 8, 1988. 
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Figure A. 1.2: RMS frequency error versus CNR for a 100-g/s trajectory 





Notice that for T a — 2 ms and Bl = 5.2 Hz, the variance is high enough at 16 dB-Hz so the 
quasi-linear approximation applies in the region where CNR > 20 dB-Hz. At point B, 50 frequency 
estimate samples are averaged to produce a frequency estimation error with variance 


E(B ! ) = ( 

n 2^r AT 
' 1=1 1=1 


• — a J — A 

1 \ 2 r At-l 1 

= ( 

-) JV^(0) + 2 2(^-|*-l|)^) 

1 ' L k=l 


(A.4) 


where N = 50 in our particular application, and RA(k) is the noise autocorrelation at point A. 
This is simply 


RA(k) = RAu>T s (kT s ) = ^ 1 [H(z)SN(z)] z=eJUlTs (A.5) 

where H(z) is the closed-loop transfer function, and Spj{z) is the power spectral density of the 
equivalent noise process (see [7] for more detail). Some typical plots are shown in Figure A.2.1.1 
for a type- III loop filter. The power spectral densities were obtained numerically and the theoretical 
autocorrelation function was obtained using an inverse FFT. This theoretical autocorrelation is very 
well approximated by an exponential function of the form 


RAu,T,(kT s ) ~ RAuT,(0)e~ 2BL ^ T ‘ 


(A.6) 


Then, Eq. (A.4) is accurately approximated by 


E(B 2 )~jR a ( 0) + 

+ 


Ua( 0) N fe-™*r- 

JV *=1 

^R A (0)J^ke~ 2B ^ 

k= 1 


(A.7) 


where i2^(0) is given in Eq. (A.6). Finally, for T s = 2 ms, N = 50, Bl = 5.2 Hz, r = 2, and k = 
0.5 


E(B 2 ) ~ 0.714^(0) (A. 8) 

Following the reasoning of [4], the variance of the frequency estimation error for the FFT is ap- 
proximated by 
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1 


(A.9) 


is 


- qou + (l - q)acR 

= +a-«) 


12 


2x 2 CAR T s 3 N(N 2 - 1) 


where q is the probability of an outlier, cr u 2 is the variance when an outlier occurs, and ocr 2 is the 
variance given by the Cramer-Rao bound. The computation of the outlier probability is identical to 
the error rate computation for the noncoherent reception of M orthogonal signals. The equivalent 
bit SNR is simply given by 


= CNR T se 
N 0 


(A. 10) 


where T ae = effective sampling interval at the input to the FFT (4 ms in our case). The equivalent 
sampling frequency is 


/« = -- = 250 Hz (A. 11) 

Ise 

The equivalent number of orthogonal signals is equal to the number of independent samples, in this 
case 25. 

For a given CNR , T se , and M, the probability of incorrent reception (of an outlier, in our terminol- 
ogy) has been tabulated in [4]. Based on this reference, and using Lagrange interpolation formulas 
for our parameters of interest, the outlier probability is easily computed. These results will be 
described in the next section. 

During operation, the FFT estimator corrects the frequency error of the CPAFC loop, so the rms 
frequency error of the FFT- CPAFC is much lower than that of the CPAFC alone. Hence, the 
variance of the estimation error at point D is essentially that of point C. 

A summary of results is shown in Figure A. 2. 1.2, which shows predicted values based on theory 
and simulation. Notice that for the system parameters previously described, and if CNR < 25 
dB-Hz, the performance of the FFT-CPAFC is actually worse than that of the CPAFC alone. This 
is simply due to the fact that for low CNRs, the probability that an outlier occurs is very high, 
producing on the average a very large estimation error. 


A. 2. 2 Dynamic Case 

When the CNR is low and there are dynamics, performance of the FFT-CPAFC loop is dominated 
by noise; therefore, the system is limited by the occurrences of outliers. The combined estimation 
error for low CNRs must then be similar to that considered in the previous section. 
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FREQUENCY, Hz 

Figure A.2.1.1: Noise spectral density 


When the CNR is high and dynamics are present, performance is dominated by the resolution of 
the FFT, suggesting a floor level. If this is the case, a simple rule of thumb to compute the variance 
of the estimation error is: 


2 ~ ( FFT bin width \ 2 _ / 250 \ 2 / 1 \ 

~ \2 • interpolation points ) D \2-32*3/ \4/ 

where Tp is the duty cycle of the dynamics, i.e., the ratio of the duration of the dynamic disturbance 
and the total observation time interval. Notice that Eq. (A. 12) is independent of the severity of the 
trajectory (provided of course that the CPAFC is still in lock, and that the input trajectory produces 
in the CPAFC loop an error in magnitude greater than the resolution of the FFT estimator). As 
an indication of how good this rule of thumb is, Eq. (A. 12) predicts a floor rms frequency error 
level of 0.65 Hz, while the simulation results predict a value of 0.40 Hz. 
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Appendix B 


Calibration Measurements 


In order to make meaningful comparisons between theoretical, experimental, and simulated perfor- 
mance, accurate calibration of the test equipment is essential. Specifically, the CNR presented to 
each algorithm must be known. Although the CNR at the output of the signal generator assembly 
(SGA) is beleived to be known to within 0.1 dB, the various losses due to cabling, attenuation, 
mixing, filtering, and A/D conversion are not well understood. Thus, the CNR of the sampled 
sequence presented to the algorithms may differ markedly from the nominal settings on the SGA. 
Therefore, we seek to develop a procedure for determining the CNR of the digital sequence directly 
from the complex samples generated by the A/D converter. Our approach will be to estimate 
the signal amplitude and noise variance first, and subsequently obtain a CNR estimate from these 
quantities. 


B.l Mathematical Model 

When operating in the PLL configuration, the GPS testbed can be modeled as in Figure B.l.l. 
The in-phase samples X{ and quadrature samples y, can be represented as 


Xi = Acos((f>i) + n Xi 


(B.la) 


Vi = Asin(&) + n yi 


(B.lb) 


where <f>i is the phase error, A the signal amplitude, and n Xi and n yi are independent zero mean 
Gaussian random variables with variance a 2, = Nq/ 2T s . Since the in-phase samples are proportional 
to A for small it is convenient to use them for estimating A and No, while adjusting the loop 
bandwidth to effectively eliminate the effects of the phase error fa. In this way, sample statistics 
can be used to estimate the signal amplitude and noise variance to virtually any desired accuracy. 
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♦i 

Figure B.1.1: Phase-locked loop block diagram 


Using the small-angle approximation cos(</>) ~ 1 - <f> 2 / 2, the mean of a;,- is 




Xi = Acos (<&) ~ A ( 1 - I = A I 1 1- 


(B.2a) 


and its variance can be obtained by means of the approximation 


cos 2 (<j>i) - cos 


if the <pi are zero-mean normal random variables with variance <r\. Assuming that the phase error 
samples are independent of the noise samples, it follows that 


var(ij) = A 2 var [cos(<fo)] + var(n X| ) ~ + a 2 (B.2b) 

For a linear PLL with loop bandwidth Bl, it is well known that if the loop bandwidth greatly 
exceeds the bandwidth of the transceiver phase process, then 


.2 _ A 0 i?£, _ Bl 
+ ~ A 2 ~ CNR 


(B.3) 


Thus, both the mean and the variance of the observed sequence depend on the variance of the 
phase error, which in turn is a function of the loop bandwidth and the carrier-to-noise ratio. We 
observe from Eq. (B.2a) that the mean value of the observed sequence can be used to estimate 
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the signal amplitude directly, provided the phase-error variance is reduced to an acceptable level. 
In principle, this can be achieved at any CNR by making the loop bandwidth sufficiently small. 
Suppose we wish to limit the bias term a 2 ^/ 2 to a value no greater than 6. Then for any CNR, the 
loop bandwidth that achieves this limit is 


B l * = 26 CNR (B.4) 

For example, at CNR = 30 dB-Hz, Bl < 2 Hz guarantees that 6 < 0.1%, while at CNR = 20 
dB-Hz, B l < 0.2 Hz guarantees the same percentage bias. As long as this bias is substantially 
smaller than the rms estimation error, its effect on the amplitude estimate can be ignored. 

The expression for the variance of the observables in Eq. (B.2b) can also be simplified when operat- 
ing with narrow bandwidths at high CNRs. Using Eq. (B.3) in Eq. (B.2b), and defining B s = l/T 5 , 
it follows that 


, x {NqBl\ f B l \ , N 0 B S 2 
var(x,) = (-j-J + — - ° (B-5) 

for Bi <C B s and Bi/CNR <1. As long as the above conditions are satisfied, we can use the 
approximations 


x, ~ A 

(B.6a) 

var(x,) ~ a 2 

(B.6b) 


These approximations help to simplify the determination of estimator performance, and will be 
assumed in the following analysis. 


B.2 Estimator Structure 


The following estimator was chosen for determining CNR from the in-phase samples: the sample 
mean is used to estimate the signal amplitude A, while the sample variance is used to determine 
No . The ratio of the mean- value squared to noise spectral level is used to estimate CNR. Thus we 
have 


N 


a= n^ x - 


t'=l 


(B.7a) 
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(B.7b) 


No = 2I > 2 = 


2T a 

N - 1 


S(*f- A) 


p 

CNR = — (B.7c) 

N 0 V ' 


If the loop bandwidth has been selected to effectively eliminate the bias term, the variance of the 
resulting estimates may be obtained by assuming that the rms estimation errors for both the signal 
amplitude and the noise spectral level are small compared to the estimates. Then we can write 


A = A + (j i (B.8a) 

No = N 0 A €n 0 (B.8b) 


CNR = 


P + 2 A€a 
No + e N 0 


~ CNR 



£Nol 

iV 0 J 


(B.8c) 


where (\ and e/v 0 are taken to be zero-mean random variables. The variance of the resulting 
estimation error may be bounded by means of the sample statistics, but care must be taken to 
account for correlations in the error terms. These correlations arise because the sample variance 
formula makes use of the sample mean. With the help of the inequality le^tyvol < 0a<?n o we can 
bound the variance of the CNR estimate as 


var - %) 5 T ^ 2 + n °~ 2(n ° 2 + lk ,° A ° N '‘ (B ' 9) 

In order to obtain numerical answers, the components of the bound must be determined. First 
consider the variance of the amplitude estimate, A. Using the model of Eq. (B.6), it follows that 
the variance of the sample mean is 


= va rM - 


(B-10) 


The variance of the sample variance involves the fourth central moment, /x 4 , of the population 
distribution. Since for zero-mean normal random variables = 3a 4 , it follows that 


= var(e^ 0 ) = 

No 


4 T 3 2 

N 



( 


N - 3 
N - 1 





8 T s 2 

N 


(B.ll) 
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for the case N > 3. Direct substitution into Eq. (B.9) yields 


var(CNR) < — [l + ( T a CNR ) 


_1 + 2 (T a CNR)~i 


(B.12) 


Referring to Eq. (B.8c), we note that the ratio CNR/ CNR = 1 + 6, where 6 is a zero-mean random 
variable whose variance is bounded by the right-hand-side of Eq. (B.12). Denoting the upper bound 
on the standard deviation of this random variable as c t u , we can approximately bound CNR/ CNR 
in the interval (1 — <r e u , 1 + <7 e “). Since o t u is itself a function of CNR, the width of these intervals 
varies as shown in Figure B.2.1, where the upper and lower envelopes are expressed in dB to 
facilitate comparison with the numerical results that follow. 


B.3 Numerical Results 


Numerical calibration results were obtained using the ARX in the PLL configuration. A loop 
bandwidth of Bi = 0.1 Hz was selected to ensure that the bias term remains less than 0.1 percent 
in the range of interest (20 < CNR < 30 dB-Hz). Since the resulting bias in the CNR estimates is 
insignificant compared to the roughly 0.2-dB estimation error envelope, the effects of bias on the 
CNR estimates can indeed be ignored. 

The calibration curve of Figure B.3.1 was obtained by processing N = 4000 in-phase samples in 
the manner described above, at 1-dB intervals in the range 20 < CNR < 30 dB-Hz. This particular 
value of N was chosen because it yielded adequate performance, and because the framework already 
existed for recording 4000 in-phase and quadrature samples (there are exactly 4000 samples in the 
8-s common trajectory used for evaluating the algorithms). Note that the spread of the data- 
points around the best-fit line is consistent with the error envelopes shown in Figure B.3.1. It can 
be concluded from the correction curve that a positive correction should be applied at all CNRs 
within the range of interest, from roughly 0.1 dB at 30 dB-Hz to nearly 0.8 dB at 20 dB-Hz. The 
reason for this miscalibration is not well understood at this time. 
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Figure B.2.1: 
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Figure B.3.1: Calibration curve 



Appendix C 


ARX Software Modifications 


This appendix describes the modifications to the ARX software that were implemented in order to 
test the various algorithms. 

In ADRCF, the subcarrier and symbol synchronization loop filters have been removed. A trajectory 
was added to the NCOl command (see Figure 5.2. 1.1). The trajectory insertion is requested by the 
operator. The operator also requests the specific algorithm for generating the NC02 command. 
The NC02 command is the frequency estimate from the algorithm chosen. These added algorithms 
are the DPLL, which is the same as that used in the carrier loop of the ARX, the CPAFC loop, 
the EKF, and the FFT-CPAFC. The ARX data acquisition was totally replaced by the GPS data 
acquisition. For the purpose of analysis, the ARX acquires I, Q, power, and SNR samples for 
different sampling intervals. These samples are continuously acquired during tracking. A tracking 
period can be any length of time. The ARX attempts to compute statistics and display them in 
real time. When the ARX computations fail to keep up with the sampling, sample acquisition is 
halted temporarily and the sample loss is noted. 

For the purpose of analysis, the data that the GPS acquires depend on the algorithm that has been 
chosen by the operator. The choice of the DPLL and CPAFC algorithms will cause acquisition of 
I, Q, NCOl command, NC02 command, the output of the cross-product, and the output of the 
cross-product integrator. The choice of the EKF algorithm will cause acquisition of I, Q, the NC02 
command, and the state vector of the Kalman filter. The choice of the FFT-CPAFC algorithm will 
cause acquisition of the NCOl command and the NC02 command. The GPS sampling interval is 
fixed at 2 ms and the period of acquisition during the track is 8 s. By assigning a 4000-element 
buffer (500 samples per second times eight seconds) for each test, no samples are lost. 

The acquisition of I and Q samples for performing the FFTs was also modified. The ARX needs 
the I and Q samples to compute FFT for preacquisition, which does not have the time constraint 
imposed by the GPS. It does, however, require that the number of FFTs performed be open ended 
so that this number can be chosen by the operator; the system always performs 80 FFTs for each 
test run. The ARX toggles between two buffers, each with the length needed to compute one FFT. 
It waits until a buffer has been freed by the program that performs the FFT before refilling that 
buffer. The GPS has 80 buffers assigned, and the last sample of any buffer is contiguous with the 
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first sample of the next buffer. The operator defines the FFT size in both systems. Additionally, 
the GPS allows the operator to set the number of data points to be transformed. 

In ARPRINT, the monitor display sections and the logging sections associated with the subcarrier 
and symbol synchronization loops were excised. Added to the monitor display were frequency 
estimation errors from the FFT and the EKF algorithm, and indicators showing which algorithm 
is currently in effect. Logic was added to control program flow associated with FFT requests and 
program flow associated with multiple test runs. A test run is one 8-s track. The GPS has five 
sections for data reduction and disk logging, as described below. 

Section 1 is for the FFT-CPAFC loop. It has two subsections depending on whether the operator 
has requested all 80 estimates of each test run or statistics. The statistics are the rms and mean 
frequency errors and lock condition for each test run and the average rms, mean, and probability 
of loss-of-lock for multiple test runs. The number of test runs are requested by the operator. The 
remaining sections are for all other algorithms. Section 2 handles the case with no dynamics. It 
generates 4000 entries of raw data and does not do multiple test runs. Section 3 handles the case J 
with no dynamics. It computes statistics and can perform multiple test runs. Section 4 handles 
the case with dynamics and generates 4000 entries of raw data. It does not do multiple test runs. 
Section 5 handles the case with dynamics. It computes statistics and can do multiple test runs. In 
Sections 2, 3, 4, and 5 there are some differences in the data reduction and disk logging capability 
depending on the algorithm in effect. 

Section 1, the output of the FFT-CPAFC algorithm, comprises 5 columns of 80 entries per column. 
Column 1 is the sample number (number 1 to 80). Column 2 is the average of 50 samples of the 
command to NC02. Column 3 is the estimate from the FFT. Column 4 is the difference between 
Columns 2 and 3. Column 5 is the sum of Columns 4 and the average of 50 commands to NCOl. 
The statistics are performed on the command to NC02 plus the command to NCOl minus the 
estimate from the FFT. 

Section 2 produces one output for the EKF algorithm and another for the CPAFC and DPLL 
algorithms. Both outputs have 4000 entries. The EKF algorithm has 7 columns and the CPAFC 
and DPLL algorithms have 6 columns. For the EKF algorithm, Column 1 is the sample number 
k for k = 1 , . . ., 4000, Column 2 is fi(k), Column 3 is xi(k), Column 4 is Column 5 is 

x$(k), Column 6 is 14 (h) and Column 7 is g(k). For the CPAFC and DPLL algorithms, Column 
1 is the sample number, Column 2 is the discriminator output, Column 3 is the integrator output 
(Columns 2 and 3 are zero for the PLL algorithm), Column 4 is the command to NC02, Column 
5 is I, and Column 6 is Q. 

Section 3 produces 3 columns of 4000 entries. Column 1 is the sample number, Column 2 is the 
command to NCOl, and Column 3 is the command to NC02. If the EKF algorithm is in effect, 
Column 2 is the average of two consecutive commands to NCOl. 

Section 4 computes and produces the statistics on both frequency estimate (the command to NC02) 
and phase estimate (Q). 

Section 5 computes and produces statistics on the error in frequency estimation. This error is the 
difference between the command to NCOl and the command to NC02. If the EKF algorithm is 
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in effect, the average of two consecutive NCOl commands is used to compute the difference. If the 
difference exceeds 250 Hz or is less than -250 Hz, an out-of-lock condition is recorded. 

In ARFAST, the averaging of FFTs that the ARX requires was removed. The 200-ms delay for 
each AP-4 access was also removed. The ARX requires a greater update rate from the A/D 
convertors than the GPS and thereby needs to reduce traffic on the multibus. The GPS requires a 
greater speed in computing FFTs. Several AP-4 commands were added so that the GPS can do a 
vector multiplication. Buffer handling was changed to conform with the changes in ADRCF. Logic 
was added to accumulate 80 FFT estimates. The power spectrum, peak search, and quadratic 
interpolation to improve the estimates are the same in both systems. 

In ARPCALC, the changes were minimal. The calculations associated with the subcarrier and 
symbol synchronization loops were simply ignored rather than removed. One gain was added for 
use in the EKF algorithm. 

In ARINITL, the initializing of the parameters used by excised portions of the software were ignored 
rather than removed. The parameters used by the real-time software are initialized and the trajec- 
tory generated. The parameter values that ARINITL reads from a disk file, upon operator request, 
were modified to reflect the added parameters of the GPS system (compatible with ARPENTR). 

In ARPENTR, the ARX portions associated with the subcarrier and symbol synchronization loops 
were excised. A section was added to accept operator inputs for algorithm selection. Should the 
operator select the FFT-CPAFC algorithm, ARPENTR asks the operator to request the number of 
I and Q samples, and the FFT size. Should the operator select the EKF algorithm, ARPENTR asks 
the operator to input the Kalman filter steady-state gains. The parameter values that ARPENTR 
writes to a disk file, upon operator request, were modified to reflect the added parameters of the 
GPS system (compatible with ARINITL). 

In ARMENUS, the statements associated with the subcarrier and symbol synchronization loops 
were excised and the statement concerning preacquisition was removed. Statements were added to 
accommodate algorithm selection. 

In ARGLOBS, the global variables not applicable to the GPS were ignored rather than removed. 
All newly defined global variables were added. 
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Appendix D 


Glossary 


s 

- second 

ms 

- millisecond 

Hz 

- Hertz 

m/s 

- meter per second 

m 

- meter 

g 

- gravity acceleration factor 

AFC 

- Automatic Frequency Control (loop) 

ALS 

- Adaptive Least Squares 

ARX 

- (NASA’s) Advanced Receiver 

BPSK 

- Bi-Phase Shift Keying 

CPAFC 

- Cross-Product Automatic Frequency Control (loop) 

CNR 

- Carrier (power) to Noise Ratio 

CRT 

- Cathode Ray Tube (terminal) 

DFT 

- Discrete Fourier Transform 

DPLL 

- Digital PLL 

DSN 

- (NASA’s) Deep Space Network 

EKF 

- Extended Kalman Filter 

FEKF 

- Frequency Extended Kalman Filter 

FFT 

- Fast Fourier Transform 

GPS 

- Global Positioning System 

IEC 

- Interstate Electronics Co. 

JPL 

- Jet Propulsion Laboratory 

MLE 

- Maximum Likelihood Estimator 

NASA 

- National Aeronautics and Space Administration 

NCO 

- Numerically Controlled Oscillator 

ODAFC 

- Overlapping DFT-based AFC 

PLL 

- Phase-Locked Loop 

RAJPO 

- (GPS) Range Applications Joint Program Office 

RMS 

- Root Mean Square 


81 



SGA 

TPS 

TTFF 

VLSI 


- Signal Generator Assembly 

- Translator- Transdigitizer Processing System 

- Time To First Fix 

- Very Large Scale Integration 
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